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PROPELLER  CAVITATION  EFFECTS  ON  THE  VESSEL  RADIATED  ACOUSTIC  SPECTRA 


I.  INTRODUCTION. 

Recently,  acousticians have  tried  to  discover  a  feasible  tech¬ 
nique  to  measure  the  prominent  source  of  acoustic  power  emanating  from  a 
merchant  vessel.  It  is  believed  that  propeller  cavitation  is  the  major 
contributor  to  merchant  vessel  radiated  noise,  but  the  fraction  of  the 
total  power  which  actually  comes  from  this  cavitation  is  not  presently 
known.  The  objective  of  this  research  effort  is  to  measure  this  percen¬ 
tage. 

Many  new  mathematical  algorithyms  are  available  to  the  experimenter 
with  the  recent  development  of  real  time  spectrum  analyzers  which  solve 
many  of  the  inherent  problems  previously  encountered  in  underwater  data 
gathering  and  processing.  Using  these  advanced  signal  processing  devi¬ 
ces,  the  coherence  function  appears  to  be  the  proper  mathematical  relatiou 
to  measure  the  effect  propeller  cavitation  has  on  the  total  acoustic 
signature. 

Once  the  mathematical  and  signal  processing  techniques  had  been 
established,  a  feasible  experimental  design  was  developed.  The  experimen¬ 
tal  setup  was  designed  with  easily  accessible  resources  to  make  the  pro¬ 
cess  repeatable,  and  also  to  allow  for  relatively  unlimited  static  data 
gathering.  Vessel  data  was  taken  at  two  locations:  in  the  near  field  of 
the  propellers,  and  in  the  far  field  of  the  entire  vessel  (see  figure  1). 
The  near  field  signal  is  totally  dominated  by  propeller  noise,  and  the  far 
field  signal  is  equidistant  from  all  noise  sources  of  the  vessel.  The 
vessel  was  an  EO  foot  displacement  craft  having  four  diesel  engines,  a  bU 
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FIGURE  1 


Hz  generator,  and  twin  three-foot  diameter  propellers.  Having  these 
characteristics  as  well  as  propellers  which  are  known  to  cavitate  at  Low 
speeds,  this  craft  is  assumed  to  be  acoustically  similar  to  a  merchant 
vessel.  The  normal  use  of  this  boat  is  to  train  midshipmen  in  seamanship 
skills  and  has  the  designation  of  Yard  Patrol  Craft  (YP).  The  YP  will 
push  against  a  seawall  allowing  for  the  procurement  of  sound  readings  on  a 
cavltating,  immobile  vessel.  Thus,  with  such  an  experimental  layout,, 
almost  unlimited  data  may  be  gathered  and  recorded  for  later  signal  pro¬ 
cessing. 

A  theoretical  investigation  was  initiated  to  show  that  the  experiment 
is  relevant  to  vessels  in  the  ocean  where  the  depth  is  essentially  infinite 
An  analysis  was  undertaken  to  determine  whether  a  mathematical  function 
could  be  used  to  essentially  eliminate  the  multi-path  and  time  delay 
problems  inherent  with  a  shallow  water  environment.  Looking  ahead,  the 
coherence  function  proves  to  be  the  ideal  algorithm  to  determine  the 
effect  propeller  cavitation  has  on  the  broadband  merchant  vessel  radiated 
acoustic  spectra. 

The  procurement  of  acoustic  data  such  as  that  described  above  is  of 
critical  Importance  to  the  Navy  today.  Presently,  underwater  acoustic 
thresholds  are  limited  by  the  amount  of  the  ambient  noise  intensity  pre¬ 
sent  at  the  hydrophone  or  detection  site.  Merchant  vessel  traffic  is 
known  to  be  the  major  source  of  ambient  noise  in  the  ocean  today. 

Therefore,  the  more  information  gained  on  the  merchant  vessel  acoustic 
spectra,  the  better  ambient  noise  is  understood,  and  hopefully  the  Navy  can 
increase  its  acoustic  thresholds  and  detect  targets  from  greater  distan¬ 
ces.  Because  of  the  major  threat  presented  by  hostile  submarines,  under¬ 
water  acoustic  research  for  the  Navy  will  always  be  given  a  high  priority. 


II.  THE  COHERENCE  FUNCTION. 

Ia  order  Co  process  Che  data  in  a  precise  and  expedient  manner,  a 
mathematical  algorithm  had  to  be  found  or  developed  which  can  be  handled 
digitally  and  also  yields  the  desired  output.  Furthermore,  the  mathema¬ 
tics  must  be  able  to  account  for  the  many  unique  qualities  involved  with 
the  acoustic  data.  The  important  qualities  include  the  multiple  paths  of 
a  water  environment  and  the  finite  propagation  time  of  sound  traveling  a 
distance  in  the  water.  Herein  lies  the  major  complexity  involved  with 
underwater  acoustic  research,  molding  the  mathematics  to  the  unique 
requirements  of  the  research.  After  investigating  many  mathematical  func 
Cions  it  was  determined  that  the  coherence  function  satisfies  all  the 
requirements  necessary  to  yield  the  appropriate  results. 

It  can  be  shown  that  the  coherence  function  measures  the  fraction  of 
the  total  power  in  the  far  field  signal  which  was  caused  by  power  present 
in  the  near  field  signal.  In  the  case  of  the  YP,  the  near  field  signal 
will  be  attained-  from  a  hydrophone  close  to  the  YP.  It  is  necessary  for 
the  near  field  hydrophone  to  be  very  near  the  propeller  because  this 
entire  signal  must  be  assumed  to  be  propeller  dominated  for  the  purposes 
of  this  analysis.  The  coherence  function  is  defined  as: 

2  I  <  Gxy  >  I  2 

r  - _ 

<'Gxx^  <Gyy^ 

where  Gxy  is  the  cross  p  wer  spe  urn,  and  G**  and  Gyy  are  simply  power 
spectral  densities.  The  power  is  defined  as: 

Gxx  -  X(f)X*(f) 


for  a  given  signal  x(t),  and  where  X(f)  is  the  Fourier  transform  of  x(t) 


In  the  experimental  set-up: 

Near  Field  Signal:  x(t)*propeller  noise 

Far  Field  Signal  :  y(t)-ax(t)+n(t)-propeller  noise+  other  noise 

where  aa  propagation  constant 


Applying  the  Fourier  transform: 

*-$x 

x(t)  *  XCfJ-Rjje 

*-$x  i(£n 

y(t)  ♦  YCfJ-aRxe  +  B^e 


Proceeding  to  the  coherence  function: 

Gxy-X(f)Y*(f)-(Rxei*X)(aRxe  **+  R^  *V 

i(<*>x“4>n> 

Gxy-aRxRx*+  «xRn*e 
<Gxy>-aRxRx 

In  the  averaging  process  inherent  in  the  coherence  function  all 
uncorrelated  terms  go  to  zero.  The  detailed  reasoning  behind  this  reduc 
tion  is  covered  in  the  time  delay  section. 

Going  on  to  the  denominator  of  the  coherence  function: 

Gxx-X(f)X*(f)-<Rxe  **)(Rxe  4>X)* 

<GXX^"RXRX 

i$n  i4>x  ^n  . 

Gyy "Y( f ) Y* ( f ) ■( aRxe  +  I^e  )(aRxe  +  R„e  ) 

<Gyy>*a  RxRx  +  ^nRn 


again  Che  uncorrelated  terms  go  to  zero  when  averaged. 
Finally, 


I  <e*y>  I  2 

Y  “  _ 

^^XX^  <tiyy> 

I  aRXRX  I  2 

2 

Y  *  _ 

<*xRx  )(a  RxRx  +RnRn  ) 

2  a  RxRx 

Y  *  _ 

a2  RxRx*+RnRn* 


y2-  INPUT  POWER 
OUTPUT  POWER 


where  the  output  power  consists  of  input  power  and  noise  power. 

Therefore,  the  above  ratio  is  exactly  what  is  needed  to  discover  just  how 
much  of  the  far  field  power  is  caused  by  propeller  cavitation. 

All  the  mathematical  operations  necessary  for  application  of  the 
coherence  function  are  performed  digitally  using  state-of-the-art  signal 
processing  techniques.  The  recorded  time  series  are  digitally  sampled  and 
separated  into  time  windows  such  that  the  Nyquist  Theorem  is  satisfied. 

The  Nyquist  Theorem  states  that  sampling  must  be  performed  at  twice  the 
highest  frequency  which  is  contained  in  the  signal.^  These  discrete 
samples  are  then  changed  from  the  time  domain  to  the  frequency  domain  via 
a  Fast  Fourier  transform  algorithm(FFT) .  This  Fourier  transform  technique 
yields  512  discrete  frequency  values  across  the  desired  frequency  range 
(0-500  Hz).  Finally,  these  frequency  values  are  multiplied,  conjugated 
and  divided  as  required  by  the  coherence  function  to  yield  the  proper 
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result.  Thi*  r-roTM*  is  carried  out  for  every  time  window.  The  duration 
of  Mrn  viMtw  La  approximately  0.8  seconds.  Consequently,  the  speed  of 
computation  ia  of  utmost  importance  in  order  not  to  lose  any  information 
from  the  incoming  time  series. 

III.  THE  EFFECT  OF  TIME  DELAY  ON  THE  COHERENCE  FUNCTION. 

From  Figure  1  it  can  be  seen  that  the  signal  received  at  the  near 
field  will  not  correspond  exactly  to  the  signal  received  at  the  same  time 
in  the  far  field.  Simply  stated,  it  takes  a  finite  length  of  time  for  the 
signal  in  the  near  field  to  propagate  to  the  far  field  in  the  underwater 
environment.  The.  question  to  be  dealt  with  is  whether  this  time  delay 
will  limit  the  effectiveness  of  the  coherence  function.  Incorporating  the 
time  delay  into  the  signal  looks  like  this: 

Near  field  Far  field 

x(t)  y(t)»ax(t-x)+n(t) 

n(t)»  noise 

Time  delay,  x,  is  the  time  it  takes  the  signal,  x(t),  to  travel 
through  the  water  to  the  far  field.  In  this  specific  case  assuming  a 
speed  of  propagation  of  1500  m/sec, 

100  meters 

T“  _  *  0.0667  seconds 

1500  meters/serond 

Applying  the  Fourier  transform  to  the  signals: 
l$x 

x(t)  ♦  Xlf)»  I^e 

l$X-iWT  4n 

+  Kne 


y(c)  ♦  YCO-aR*. 


Proceeding  to  the  definition  of  the  coherence  function  once  again,  the 
numerator  reduces  to: 


iy  X  i$  x  “iWT  q 

Gxy-X(f)Y*(f)-(Rxe  KaRxe  e  +  ^e  )* 

io>T  i(<t>x'<f>n) 

GXy*aRxR-x*e  +  R^R^e 


Considering  each  term  separately  as  vectors  in  phase  space: 
aRxRx*eimi  -  always  has  the  same  phase,  tux,  so  when  averaged  the  vectors 
tend  to  add  coherently. 


f  A  coherent  vector  sum. 


*  iC'i’x-^n) 

RxRn  e  -  since  x  and  a  are  independent  processes,  the  phases  ot 

the  vectors  vary  randomly  over  the  average,  and  hence 
the  term  is  insignificant  when  compared  to  the  first 


A  random  vector  sum. 

Herein  lies  the  explanation  of  why  the  averaging  process  within  the 
coherence  function  is  so  helpful  in  this  application.  Only  the  vectors 
which  sum  coherently  are  significant  in  the  coherence  measure,  and  these 


same  vectors  represent  the  coherent  power  in  the  far  field. 


The  denooinacor  of  Che  coherence  is  unaffected  by  Che  time  delay 


because  Che  cross  power  spectrum  of  che  same  signal,  or  Gyy,  negates 
the  phase  difference.  In  the  average,  the  vectors  representing  both  Gxx 
and  Gyy  have  the  same  relative  phase  difference,  and  therefore  tend  to 
sum  coherently. 

<  Gxx>*<x<  £  y x*  ( ■ f )  >■  i &xRx* 

<Gyy>-<Y(f)Y*(f)>-a2RxRx*+  R^* 


Cross  terms  are  incoherent  and  go  to  zero  when  averaged  as  described  earlier 

<GxxXGyy>-(RxRx*)Ca2  EhtHx*+RnRn*) 

Therefore,  che  coherence  now  looks  like: 

y2m  |  aRxRx*eiwT  I  2  „  a2jtxRx* 

(  RXRX*  )  (  a2  RxRx*+RnRn* >  *  RxRx*+  RnRn* 

Y2-  INPUT  POWER _  -  INPUT  POWER 

INPUT  POWER  +  NOISE  POWER  OUTPUT  POWER 

Once  again  che  desired  proportion  of  che  output  power  caused  by  input 
power  in  che  far  field  is  obtained.  Consequently,  che  time  delay  does  not 
limit  Che  effectiveness  of  the  digitally  processed  coherence  function. 

IV.  DIGITAL  TIME  WINDOWS. 

A  major  problem  involved  in  the  digital  processing  of  the  time 
delayed  data  is  Che  adjustment  of  the  sampling  windows  to  get  che  desired 
"overlap'*  of  che  two  time  series.  The  reason  this  is  important  lies  in 
che  fact  chat  both  che  near  and  far  field  signals  are  recorded  simulta- 
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neously  on  an  analog  cape  recorder.  When  Che  signal  processor  diglcally 
samples  Chls  recording  ic  performs  all  ics  corapuCaCions  with  chese 
parallel  windows.  In  Chls  experiment,  a  faulty  coherence  measure  would  be 
gained  because  Che  far  field  signal  of  inceresc  does  noc  run  contiguously 
wich  Che  corresponding  near  field  signal.  Even  chough  Che  propagation 
time  is  so  small,  some  sort  of  digital  delay  must  be  implemented  to  obtain 
an  exact  coherence  measure.  In  Che  frequency  band  of  inceresc  for  this 
experiment,  0-500  Hz,  Che  propagation  Cime  of  approximately  0.07  seconds 
accounts  for  only  10  percent  of  the  0.8  second  digital  time  window.  A 
typical  digital  time  window  display  with  insufficient  time  delay  would 
look  like  this: 


Near  field 


Far  field 


T-  time  delay  due 
to  propogatlon 


Low 

coherence 


As  mentioned  earlier,  the  entire  signal  in  the  near  field  of  the  pro¬ 
peller  is  considered  to  be  the  desired  cavitation  noise.  Therefore,  all 
of  the  far  field  time  window  which  is  not  overlapped  with  the  propeller 
noise  signal  emitted  in  the  near  field  is  considered  incoherent  noise. 

The  less  the  amount  of  overlap  the  more  inaccurate  the  measure  of  coherent 
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power  la  Che  far  field*  Fortunately,  Che  signal  processor  contains  a 
built-in  digital  cine  delay  function  which  allows  for  Che  adjustment  of 
Che  time  windows  for  maximum  overlap.  The  signal  processor  actually 
samples  Che  near  field  signal  and  Chen  wales  0.07  seconds,  the  programmed 
ciae ,  before  sampling  the  far  field  signal: 


Near  field 


V  J  I 


V 

-r-*i 


Maximum  coherence  where: 


Far  field 
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•"x 

time 

__  \ _ -  -  ■ 

delay 

time  delay-propagation 

time 


Thus,  Che  properties  of  diglcal  time  windows  allow  for  the  manipula¬ 
tion  of  Che  gathered  data  in  order  to  ensure  that  the  derived  results  are 
accurate. 
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V.  THE  EFFECT  OF  MULTI-PATH  ON  THE  COHERENCE  FUNCTION. 

A  final  theoretical  determination  must  be  undertaken  to  see  if  the 
coherence  function  is  adversely  affected  by  the  shallow  water,  multi-path 
environment  present  at  the  experimental  site.  The  true  question  being 
posed  here  is,  can  the  more  easily  used  shallow  water  environment 
correctly  represent  the  deep  water  environment  of  interest  to  sonar 
acousticians.  Plctorially  the  multi-path  situation  looks  like: 
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Where  t {  and  t±  are  the  two  propagation  times  for  the  different  paths. 

Again  using  the  Fourier  transform  (the  spectrum  uses  the  Fast  Fourier 
transform,  FFT) : 


x(t)  ♦  X(f)-Rxe 

i$x-iuxi  i$x-iux2  i$n 
y(t)  +  YCfJ^aRjje  +bRxe  +Rne 


a  and  b  are  again  propagation  factors. 

Applying  these  frequency  spec t rums  to  the  coherence  function,  the  numera¬ 
tor  becomes: 

i$x  i$x-iux  [  i$x— iux  2  l^n 

Gxy»X(  f )  Y*(  f  )■(  Rxe  )(aRxe  +bRxe  +Rne  )* 

•  • 

iojXi  h*)X2  U-frx-^n) 

Gxy”aRxRx  e  +bRxRx  e  +RxRn  e 

The  Rx^n*  term  goes  to  zero  as  the  average  is  taken,  this  vector  is  ran¬ 
dom,  whereas  the  other  vector,  KXRX*,  adds  coherently  in  the  average,  so: 

iux 1  iux  2 

<GXy>-aRxRx*e  +bRxRx*e 

Continuing  on  with  the  numerator, 

_  iux  |  iux  2  iux  1  iux  2 

|  «Jxy>  |  »  (aRxRx*e  +bRxRx*e  )(aRxRx*e  +bRxRx*e  )* 

iu(Tl-X2)  ho  (  x  2~t  1 ) 

I  <Gxy>  |  -(aRxRx  )  +(bRxRx  )  +ab(RxRx  )  +e  ) 

|  <Gxy>  I  2-(aRxRx*)Z+(bRxRx*)2+2ab(RxRx*)ros(uAx) 


where  Ax*x^-x2 


■Hi  CU 


The  denominacor  becomes: 


^Ojm>“X(f)X  (f)-RxRx 

<Gyy>«a2  RxRx*+b2  RxRx*+abRxRx*cos(uAT ) +knRn* 


Once  again  Che  cross  terms,  RjfR^,  are  neglected  in  che  averaging  process. 
Finally,  che  denominacor  becomes: 


<Gxx><Gyy>-RxRx*(a2RxRx*+b2RxRx*+abRxRx*(2c°s(a>AT  )  )+RnRn*) 


The  coherence  is: 


a2  RxRx*+b2  RxRx*+2abRxRx*cos(wAx ) 


a2  RxRx*+b2  RxRx*+2abRxRx*cos  (wAt )+RnRn* 


y2-  INPUT  POWER _  -  INPUT  POWER 

INPUT  POWER  +  NOISE  POWER  OUTPUT  POWER 


Thus,  che  multi-path  environment  does  not  represent  a  problem  in 
determining  the  coherence  value  because  che  desired  ratio  of  input  power 
to  output  power  is  again  obtained.  Although,  che  multi-path  problem  is 
much  more  prominent  in  a  shallow  water  environment,  the  ocean  also 
portrays  multi-path  characteristics.  Acoustic  power  may  propagate  via 
many  differing  paths  such  as  convergence  zones  and  the  shallow  waveguides 
above  the  thermal  layer. 


VI.  SUMMARY  OF  THEORETICAL  ANALYSIS. 

After  a  very  general  mathematical  analysis  of  the  problems  associated 
with  digitally  applying  data  from  a  time  delayed,  multi-path  environment, 


the  coherence  function  appears  to  be  the  most  robust  operation  available 
for  obtaining  the  desired  ratio.  Two  unique  qualities  inherent  in  the 
coherence  function  which  make  it  very  useful  in  this  specific  application 
are  the  averaging  process  and  the  frequency  dependence.  The  averages 
taken  within  the  function  basically  integrate  out  the  cross  terms  which 
are  uncorrelated.  The  spectrum  analyzer  also  can  be  adjusted  to  take  dif¬ 
ferent  numbers  of  averages  to  gain  a  quicker  (fewer  averages)  or  more 
accurate  (more  averages)  result  depending  on  the  requirements  of  the 
experiment.  When  more  averages  are  taken  the  phase  random  vectors  become 
more  insignificant  and  the  coherent  vectors  are  overwhelming  dominant. 

The  high  degree  of  precision  desired  necessitates  the  long  period  of  data 
gathering  which  is  incorporated  into  the  experimental  design. 

The  narrow-band  frequency  dependent  characteristics  of  the  coherence 
function  are  actually  the  underlying  principles  which  allow  for  its  appli¬ 
cation  under  the  specific  requirements  of  this  experiment.  When  operating 
in  the  frequency  domain,  any  dominant  tonal  qualities  will  only  affect  a 
small  frequency  bin,  whereas  in  the  time  domain  any  strong  tone  would  mask 
any  truly  useful  information.  Because  of  the  large  amount  of  energy  loca¬ 
lized  at  60  Hz  at  the  experimental  site,  a  function  which  operates  in  the 
frequency  domain  is  a  necessity;  that  is  the  60  Hz  and  harmonics  would 
tend  to  dominate  any  broadband  measure.  A  coherence  value  is  computed 
within  the  spectrum  analyzer  for  each  of  the  512  frequency  bins.  Thus, 
any  frequency  dependence  in  the  propagation  time  or  the  multi-path  charai — 
terlstics  would  be  treated  on  a  frequency  by  frequency  basis.  Other  time- 
dependent  mathematical  functions,  such  as  the  cross  correlation  function, 
would  not  reveal  any  frequency  dependent  signal  characteristics. 
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Therefore,  Che  coherence  function  was  the  mathematical  operation  chosen  to 
describe  the  proportion  of  broadband  merchant  radiated  noise  which  is 
caused  by  propeller  cavitation. 

VII.  THE  EXPERIMENT  ITSELF. 

Data  was  gathered  as  described  in  the  introduction  for  approximately 
twenty  minutes  at  various  shaft  speeds.  The  near  field  hydrophone  was 
kept  within  five  feet  of  the  propeller  at  all  times,  while  the  far  field 
hydrophone  was  deployed  approximately  350  feet  from  the  vessel.  The  350 
foot  distance  allows  for  all  sources  of  acoustic  intensity  emanating  from 
the  SO  foot  YP  to  be  considered  as  a  point  source  of  radiation.  Ambient 
noise  level  readings  from  the  far  field  hydrophone  (figure  2)  show  that 
the  60  Hz  tone  and  its  harmonics  dominate  the  ambient  noise  in  the  Severn 
River.  This  abundance  of  60  Hz  power  is  expected  in  the  Severn  because  of 
the  many  power  and  communication  lines  along  the  bottom  of  the  river.  The 
remainder  of  the  ambient  noise  spectrum  is  "bumpy-white"  or  Gaussian,  as 
it  would  be  in  the  ocean. 

Figure  3  displays  the  signal-to-noise  ratio  created  by  the  ravitating 
YP's  acoustic  spectrum  in  the  far  field.  As  can  be  seen  from  the  figure, 
the  signal  emanating  from  the  YP  was  predominantly  over  the  ambient  noise 
in  the  120-300  Hz  range.  The  remainder  of  the  far  field  spectrum  with  the 
YP  operating  had  no  perceptible  difference  from  the  ambient  levels  in  the 
0-500  Hz  range.  During  the  data  gathering  process,  a  large  amount  of 
cavitation  was  visible  from  the  stern  of  the  YP.  To  verify  whether  this 
cavitation  was  truly  the  dominating  source  of  the  acoustic  intensity 
received  in  the  far  field,  the  coherence  function  was  applied  to  the  near 
and  far  field  signals.  Figure  4  displays  the  coherence  function  from 
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O-500  Hz,  64  time-delayed  spectral  frames  were  averaged  to  produce  an 
accurate  representation.  The  64  averages  provide  a  very  high  confidence 
limit  of  approximately  90  percent  which  simply  indicates  that  all  phase 
and  magnitude  information  will  be  accurate  to  within  10  percent.  The 
time-delayed  signal  accounts  for  the  time  for  the  signal  to  propagate  from 
the  near  field  to  the  far  field. 

As  expected,  the  coherence  values  in  the  120-300  Hz  range  were  much 
higher  when  compared  to  the  other  frequencies.  These  high  coherence 
values  verify  that  a  high  proportion  of  the  far  field  intensity  at  these 
frequencies  was  caused  by  propeller  cavitation.  Because  propeller  cavita¬ 
tion  is  a  modulated  broadband  Gaussian  phenomenon,  the  coherence  values  in 
•• 

the  region  where  good  signal-to-noise  existed  are  expected  to  be  flat 
(constant)  and  not  contain  all  the  peaks  and  valleys  as  in  figure  4.  Of 
note  is  the  fact  that  many  of  these  peaks  and  valleys  in  the  region  of 
interest,  120-300  Hz,  correspond  to  frequencies  where  there  were  either  . 
high  or  low  signal-to-noise  levels  in  the  far  field  (figure  3).  Later,  it 
will  be  discussed  how  this  "jaggedness"  in  the  coherence  cannot  be  totally 
explained  using  signal-to-noise  considerations. 

Another  interesting  feature  of  the  coherence  values  is  the  fact  that 
it  Is  very  nearly  zero  at  all  the  harmonics  of  60  Hz.  During  the 
averaging  process,  the  coherence  values  for  these  tones  decreased  uni¬ 
formly,  while  a  majority  of  the  other  frequencies  displayed  a  somewhat 
oscillatory  nature.  This  unusual  behavior  is  characteristic  of  two 
strongly  uncorrelated  sources  of  60  Hz  power.  Taking  many  factors  into 
account,  it  was  determined  that  t.ie  two  sources  of  60  Hz  acoustic  power 
are  the  ambient  intensity  which  was  already  mentioned  which  dominates  the 
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far  field  hydrophone,  and  the  60  Hz  power  emanating  from  the  generator  on 
board  the  YP.  As  expected,  these  two  sources  of  60  Hz  intensity  are  very 
much  uncorrelated  (ie.  random  phase  relationship),  and  thus  produce  this 
very  low  coherence. 

This  interesting  behavior  of  the  coherence  function  at  60  Hz,  which 
was  first  thought  to  be  an  error  in  the  processing,  actually  paved  the  way 
to  the  only  feasible  explanation  for  the  peaks  and  valleys  in  the 
coherence.  The  valleys  in  the  coherence  function  in  the  120-300  Hz  range 

can  be  considered  analogous  to  the  minimums  at  60,120,180 . 480  Hz, 

there  must  be  some  source  other  than  propeller  cavitation  which  dominates 
the  far  field  spectrum.  This  indicates  that  not  all  the  signal-to-noise 
as  seen  in  figure  3  between  120  and  300  Hz  is  caused  by  propeller  noise. 

To  verify  these  hypotheses  concerning  the  effect  of  propeller  cavita¬ 
tion  on  the  broadband  vessel  radiated  acoustic  spectra,  many  theoretical 
and  experimental  investigations  were  undertaken.  As  a  result  of  this 
first  set  of  data  it  was  determined  that  more  data  should  be  gathered  to 
discover  the  nominal  value  of  the  coherence  in  the  given  experimental 
environment.  The  puzzling  problem  was  the  fact  that  from  a  purely  mathema¬ 
tical  analysis,  a  signal-to-noise  level  of  20  dB  (i.e.  150  Hz-  figure  3) 
corresponds  to  a  coherence  of  0.99999  and  from  the  experimental  data  the 
coherence  is  only  0.9  (figure  4). 

VIII.  DETERMINING  THE  MAXIMUM  EXPECTED  COHERENCE. 

To  determine  the  nominal  value  of  the  coherence,  the  YP  was  replaced 
by  a  tonal  projector.  Tones  were  examined  which  had  high  signal-to-noise 
in  the  far  field,  and  a  coherence  of  .9  was  observed.  Thus,  0.9  seems  to 
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be  Che  maximum  coherence  supportable  by  this  acoustic  environment.  The 
same  positioning  of  the  hydrophones  was  used  but  instead  of  the  YP,  a 
transducer  was  used  as  the  near  field  source  of  acoustic  energy.  The 
transducer  was  powered  by  a  tone  generator  which  emitted  one  tone  every 
third  of  an  octave  between  20  and  500  Hz  (15  tones  total).  Knowing  the 
exact  frequencies  of  these  generated  tones  allows  for  the  determination  of 
characteristics  of  the  coherence  function  based  on  these  representative 
frequencies. 

As  seen  from  the  power  spectrums  of  the  near  field  and  the  far  field 
signals  (figure  5),  the  tones  were  very  much  attenuated  in  the  shallow 
water  environment.  For  a  few  of  the  tones,  although,  enough  signal-to- 
noise  existed  in  the  far  field  to  yield  fairly  high  coherence  values 
(figure  6).  Maximum  values  in  Che  .8-. 9  range  for  the  250  Hz  and  490  Hz 
tones  verify  that  the  values  gained  with  the  YP  were  indeed  maximums.  The 
reason. for  the  decrease  in  the  signal-to-noise  below  approximately  luO  Hz 
is  that  the  cut-off  frequency  in  the  experimental  environment  is  very  near 
this  value.  The  cut-off  frequency  is  that  frequency  where  propagation 
decreases  rapidly  on  account  of  the  waveguide  properties  of  the  surface 
and  the  boctom  of  the  river.  The  reason  more  of  the  tones  above  100  Hz 
did  not  yield  a  higher  coherence  is  the  fact  that  the  necessary  power 
could  not  be  put  into  these  tones  in  order  to  dominate  the  ambient  noise. 

A  point  of  interest  Involved  with  the  coherence  function  derived  from 
this  experiment  is  the  very  high  coherence  values  at  bO  Hz  and  its  har¬ 
monics,  i.e.  both  hydrophones  are  dominated  by  the  same  power  cable.  This 
result  is  exactly  opposite  the  results  gained  from  the  YP  data  where  the 
coherence  values  were  very  near  zero.  This  discrepancy  between  the 
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coherence  values  aC  these  frequencies  actually  led  to  the  conclusion  that 
there  oust  be  another  dominant  source  of  acoustic  power  on  board  the  YP 
other  than  propeller  noise,  particularly  sources  not  at  the  harmonics  of 
60  Hz.  With  the  tone  generator,  a  60  Hz  tone  was  not  emitted  as  it  was 
from  the  YP,  and  therefore  the  major  source  of  60  Hz  power  comes  from  the 
ambient  noise.  Since  this  same  ambient  noise  energy  is  predominant  in 
both  the  near  and  far  field  signals,  the  power  at  these  frequencies  should 
be  highly  coherent. 

Another  cross  check  of  the  coherence  function  was  performed  from  the 
data  gathered  during  this  second  experiment.  Simultaneous  ambient  noise 
spectrums  were  taken  from  both  near  and  far  field  hydrophones  and  the 
resulting  coherence  displayed  minimal  coherence  across  the  frequency  range 
(figure  7).  This  low  coherence  is  expected  because  ambient  noise  is  pre¬ 
dominantly  white  or  Gaussian  noise  which  is  always  random  and  thus  inco¬ 
herent  with  other  signal  sources.  The  60  Hz  tone  and  its  harmonics 
display  high  coherence  in  this  instance  as  with  the  data  gathered  with  the 
tone  generator  because  of  the  coherent  properties  of  these  special  fre¬ 
quencies  in  the  ambient  spectrum. 

From  this  secondary  experiment  many  useful  characteristics  of  the 
coherence  function  have  been  observed.  First  of  all,  a  nominal  value  of 
the  coherence  is  approximately  0.9,  so  that  the  value  gained  from  the  YP 
data  is  indeed  a  maximum.  Secondly,  the  60  Hz  coherence  characteristics 
suggested  the  possibility  that  propeller  cavitation  is  not  the  sole  source 
of  vessel  radiated  acoustic  energy.  Finally,  the  verification  of  the  ran¬ 
dom  properties  of  ambient  noise  were  demonstrated  when  the  two  ambient 
time  series  were  placed  into  the  coherence  function. 
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IX.  COHERENCE  FUNCTION  CHARACTERISTICS. 


From  figure  3,  1C  can  be  seen  that  Che  coherence  function  has  very 
abrupt  discontinuities  throughout  the  frequency  range.  When  comparing  the 
marked  increases  and  decreases  in  the  coherence  values  with  the  signaL-to- 
noise  levels  (figure  2),  an  apparent  correlation  was  noted  between  the 
two.  From  the  mathematical  analysis  performed  earlier,  although,  the 
sharp  changes  cannot  be  explained  within  the  parameters  established  in  the 
analysis.  Since  a  purely  signal- to-noise  approach  cannot  explain  the 
jaggedness  of  the  coherence  there  must  be  some  other  cause  for  this  beha¬ 
vior  of  the  function. 

Two  possible* explanations  for  this  characteristic  are  the  presence  of 
a  near  field  noise  source  or  possibly  the  presence  of  far  field  power 
coming  from  parts  of  the  vessel  other  than  the  propeller  on  board  the  YP. 

A  near  field  noise  source  implies  that  the  propeller  noise  is  not  comple¬ 
tely  dominating,  and  the  mosc  likely  cause  of  this  unwanted  power  is  from 
flow  noise  around  the  hydrophone.  A  slight  amount  of  water  flow  was  suc¬ 
tioned  past  the  stern  of  the  YP  into  the  propellers,  so  the  possibility  of 
flow  noise  must  not  be  neglected.  Since  the  deployed  hydrophone  displayed 
little  deflection  due  to  these  currents,  and  no  bubbling  was  observed 
about  the  hydrophone,  the  possibility  of  flow  noise  was  not  too  great.  To 
verify  that  flow  noise  actually  was  not  a  contributor  to  the  near  field 
signal,  a  mathematical  treatment  was  undertaken. 

The  second  and  more  probable  cause  for  the  jaggedness  of  the 
coherence  is  a  source  of  far  field  power  coming  from  sources  on  the  YP 
other  than  the  propellers.  If  this  were  the  case,  the  signal-to-noise 
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Levels  (figure  2)  would  not  completely  represeat  the  power  in  the  far 
field  emanating  from  the  propeller  which  was  assumed  earlier.  Some  of  the 
signal-to-noise  level  in  the  far  field  could  be  attributable  to  the  diesel 
engines  or  the  generator  on  board  the  YP.  Naturally,  power  from  these 
sources  would  not  be  coherent  with  the  near  field,  propeller  signal.  If 
indeed  sources  other  than  Che  propellers  are  dominating  the  far  field 
radiated  spectrum,  this  would  be  contrary  to  the  present  belief  that  the 
merchant  spectra  is  completely  dominated  by  propeller  cavitation  in  the 
far  field. 

The  following  mathematical  analysis  is  performed  to  determine  if 
enough  flow  noise  existed  in  the  near  field  to  account  for  discontinuities 
in  the  coherence  function.  from  the  original  analysis: 

y2 -INPUT  POWER  (in  far  field) _ 

INPUT  POWER  (in  far  field)  +  NOISE  POWER 


Sf*+N* 


From  figure  2,  the  signal-to-noise  is  20  dB  at  150  Hz,  so: 

Far  field  signal  •  20  dB  +  ambient  noise 
Sf  -  20  +  N 

The  associated  powers  are  related  by: 

Input  power  (in  far  field)  -  100  *  Ambient  noise  power 
Sf*  -  100  *  N* 


The  coherence  now  looks  like: 


Y2-  Sf*  -  10QN*  -  100  -  0.99 

-  101 

Sf*+N*  1Q0N*+N* 

Maximum  coherence! 

Previously,  a  coherence  value  of  0.9  was  determined  to  be  the  nominal 
maximum  for  the  given  experimental  environment.  At  150  Hz,  in  figure  3, 
the  coherence  is  in  fact  0.9,  or  a  maximum.  The  mysterious  feature, 
though,  is  the  0.6  drop  in  the  coherence  between  150  and  155  Hz,  from  0.9 
to  0.3.  A  signal-to-noise  approach  cannot  explain  this  abrupt  change: 

Signal-to-noise  at  155  Hz  is  5  dB,  or 

Far  field  signal  (Sf)  *  5  dB  +  Ambient  Noise  (N) 

and.  Input  power  (Sf*)  ■  3.16  *  Noise  power  (N*) 

Finally  the  coherence  value  becomes: 

Sf*  (3.16)N*  3.16 

Y  “  _  ■ _  ■ _  ■  0.76 

Sf*  +  N*  (3. 16)N*  +  N*  4.16 

Thus,  the  signal-to-noise  level  at  155  Hz  accounts  for  only  a  0.2  drop  in 
the  coherence.  The  next  task  was  to  mathematically  check  the  possibility 
of  a  flow  noise  contribution  to  the  near  field  signal: 
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Near  field 

x(t)+n(t) 

n(t)«.  flow  noise 

^♦x  *+ni 

XCf)-^  +  Rn2e 

X(f)-Sn+N1 

The  coherence  funcClon  becomes: 

Y 


Far  field 

y(t)-ax(t)  +  n(c) 

i<|>x-iwT  i(j>n 

Y(f  J-aRjje  +1^  * 

Y(f)-Sf+N2 

a2(RxRx*)Z 


a  (Rx^x  J+RxRx  (Rn2®n2*+a  *Rn2Rn2* 


and  reduces  Co: 

2  sn*  Sf* 

Y  * _ 

(Sn*+Nl*><Sf*+N2*) 

Y2-  (POWER  SIGNAL-NEAR  FIELD) (POWER  SIGNAL-FAR  FIELD) 

(TOTAL.  POWER  NEAR  FIELD) (TOTAL  POWER  FAR  FIELD) 

Fortunately,  this  equation  reduces  to  the  same  as  that  derived  with  no  flow 
noise  when  the  near  field  signal  (Sn*)  dominates  the  flow  noise  (N*)  and  the 
near  field  terms  cancel* 


Y2-  POWER  SIGNAL-FAR  FIELD 
TOTAL  POWER  FAR  FIELD 

Just  dealing  with  the  portion  of  the  coherence  function  containing  the  near 
field  terms,  and  treating  it  as  an  adjustment  factor: 


Coherence  adjustment  factor 


Sn*+Ni* 


Using  various  assumed  signal-to-noise  levels  In  the  near  field  the  following 
table  was  constructed: 


Signal-to- 
nolse  ratio 
(Sn-N!) 


Power 

ratio 

(Sn*/N,*) 


Coherence  adjustment 
factor 

(SnVCSn^Ni*) 


This  adjustment  factor  would  decrease  the  coherence  by  the  calculated 
value  If  flow  noise  did  exist  In  the  near  field.  As  verified  above,  if 
the  near  field  signal  is  more  than  10  dB  above  the  flow  noise  contribu¬ 
tion,  the  adjustment  factor  Is  essentially  unity  and  has  negligible  effect 
on  the  coherence.  Flow  noise  is  extremely  white  in  nature  and  should 
therefore  affect  all  frequencies  in  a  like  manner.  Consequently,  the 
coherence  adjustment  factor  should  reduce  all  the  coherence  values  the 
same  relative  amount,  and  thus  fails  to  explain  the  0.6  drop  between  150 


and  155  Hz.  There  is  a  possibility  that  the  flow  noise  is  highly  fre¬ 
quency  dependent,  and  enough  of  its  power  is  concentrated  at  155  Hz  to 
account  for  the  0.6  drop.  For  this  to  be  the  case,  though,  the  flow  noise 
would  have  to  completely  dominate  the  near  field  propeller  noise  at  this 


frequency  which  is  very  highly  unlikely  with  the  cavitating  propeller  only 
five  feet  away  from  the  hydrophone. 

The  discontinuous  characteristics,  therefore,  must  be  explained  by 
the  presence  of  far  field  signal  sources  other  than  the  propeller  noise. 

It  was  this  jaggedness  that  was  first  considered  an  experimental  error 
which  eventually  substantiated  the  analysis  of  the  coherence  as  the  actual 
measure  of  the  effect  propeller  cavitation  has  on  the  far  field  merchant 
radiated  acoustic  spectra.  Therefore,  the  sharp  discontinuties  in  the 
coherence  indicate  that  sources  of  acoustic  power  from  the  vessel  other 
than  the  propeller  are  dominating  the  far  field  spectra  at  certain  fre¬ 
quencies.  This  is  contrary  to  the  present  beliefs  that  the  merchant 
vessel  far  field  spectrum  is  totally  dominated  by  broad-band,  Gaussian, 
propeller  cavitation. 

X.  SUMMARY  OF  RESULTS. 

Conclusive  evidence  has  been  gained  which  shows  that  propel Ler  cavi¬ 
tation  is  not  the  dominant  source  of  acoustic  energy  emanating  from  a 
merchant  vessel  over  many  sonar  frequency  bands.  Even  with  the  YP,  pro¬ 
peller  cavitation  energy  did  not  even  dominate  a  vessel  which  had  no  flow 
noise  about  the  hull  (which  is  a  common  contribution  to  the  far  field 
acoustic  intensity).  Without  this  flow  noise,  propeller  cavitation  would 
be  expected  to  enormously  overwhelm  the  other  sources  of  noise  onboard  the 
vessel.  At  many  frequencies  where  good  signal-to-noise  existed,  a  low 
coherence  was  gained  indicating  a  non-propeller  dominated  portion  of  the 
far  field  spectrum.  These  findings  about  propeller  cavitation  energy 
which  is  propagated  great  distances  can  hopefully  begin  to  characterize 


Che  qualities  of  anbienc  noise  in  Che  ocean.  Through  this  charac¬ 
terization,  Navy  sonars  will  be  able  to  detect  enemy  targets  (submarines) 
at  ouch  greater  distances  and  with  much  greater  accuracy. 

The  data  for  the  YP  should  very  nearly  represent  data  from  a  standard 
merchant  in  the  deep  ocean.  The  many  possible  paths  of  propogation  in  the 
Severn  River  is  quite  similar  to  the  numerous  propagation  paths  in  the 
ocean  such  as  convergence  zones  and  the  shallow  waveguide  effect.  Also, 
the  YP  is  essentially  a  scaled  down  version  of  a  standard  merchant,  a 
majority  of  oceangoing  merchants  now  have  diesel  engines,  highly  cavi- 
tating  propellers  as  well  as  electrical  generators  (50  hz  or  60  Hz). 
Considering  these_  factors,  Che  data  and  the  resulting  findings  can  be 
transfered  to  the  blue-water  environment  of  interest. 

With  the  more  complex,  real-time  digital  spectrum  analyzers  being 
marketed  today,  both  the  quantity  and  quality  of  underwater  acoustic  data 
will  be  enhanced  Co  a  large  extent.  With  the  increase  in  precision  and  in 
speed  due  to  these  new  breeds  of  computers,  nearly  all  of  the  previously 
disastrous  effects  can  be  handled  with  relative  ease.  Problems  such  as 
doppler  broadening,  multi-path  interference  and  precision  hydrophone  and 
projector  placement  can  basically  be  eliminated  via  digital  processing 
techniques.  The  future  for  underwater  acoustic  research  looks  bright,  and 
as  long  as  there  are  submarines  in  the  oceans  Che  Navy  need  wii,  always 
exist.  Therefore,  this  project  answers  a  portion  of  a  problem  which  is 
of  importance  to  the  Navy  as  well  as  offering  an  opportunity  to  be  exposed 
to  an  exciting  field  of  physics. 
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